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of transmission resonances. In contrast, for a parallel double-dot we find parameter regimes in which 
the two-particle interaction generates additional resonances. 
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I. INTRODUCTION 

Electronic transport through multi-level quantum dots 
and several coupled dots is currently of great ex- 
perimental and theoretical interest due to the pos- 
sible application of such systems in interferometers 
and for charge- and spin-based quantum informa- 
tion processing! ^i^i^i'^i^i^i'^'i^i^i-^^i-^^i^^ The smallness of the 
mesoscopically confined electron droplet — the quantum 
dot — leads to fairly large energy level spacings and at suf- 
ficiently small temperatures T only a few levels must be 
considered. The physics is strongly affected by charging 
effects^ (Coulomb blockade) and local Coulomb correla- 
tions, e.g. leading to the Kondo effect .iii^J^iiSiiLia 

We here present a new method to describe the low- 
temperature transport through both a single multi-level 
dot and several coupled dots. To keep the notation short 
when referring to both situations in the following we de- 
note the dot(s) as the mesoscopic system. We investigate 
the linear response conductance G through the setup as 
a function of a gate voltage Vg that shifts the energy lev- 
els. The connection between the mesoscopic system and 
the leads (the reservoirs) is modeled by low-transmission 
tunneling barriers. 

The transport through a single-level dot with spin de- 
generacy (see Fig. ^A) is well understood. For simplic- 
ity we only consider equal couplings to the left and right 
lead. At small T and for noninteracting dot electrons 
G{Vg) shows a Lorentzian resonance of unitary height 
2e^ /h placed roughly at the energy of the single-particle 
level of the dot. The full width 2r of the resonance sets 
an energy scale F that is associated with the strength of 
the tunneling barriers. Including a Coulomb interaction 
U between the spin up and down dot electrons the line 
shape is substantially altered as can be seen from the ex- 
act T — Bethe ansatz solutionr 8 For increasing U /F it 



is gradually transformed into a box-shaped resonance of 
unitary heights with a plateau of width U and a sharp 
decrease of G to the left and right of itpiiiSiSSiSi For gate 
voltages within the plateau the dot is half-filled implying 
a local spin-1/2 degree of freedom on the dot. Thus, the 
Kondo effeciji^ leads to resonant transport throughout 
this Kondo regime. A linear chain of a few coupled single- 
level dots (see Fig. ^ B) presents a simple extension of 
the single-dot case and was studied recentlyi^SiSi^ This 
system can equally be viewed as a short Hubbard chain. 

In more complex cases in which electrons can pass the 
mesoscopic system following different paths, not only cor- 
relations but also quantum interference plays an impor- 
tant role, leading e.g. to the Fanc^^ effect. Realizations 
of such a situation that have been studied theoretically 
are: i) A single dot coupled to two leads that are in addi- 
tion coupled by a direct transmission channel (Aharonov- 
Bohm interferometer with an embedded dot) as sketched 
in Fig. n]C.^® ii) A system of two dots, one directly con- 
nected to the two leads and the other side-coupled to 
the first oneSLSSiSS shown in Fig. n]D. in) Parallel dou- 
ble dots coupled to common leads'^"'^^-'^^ '^'^ (see Fig-^E) 
that have also been realized experimentallyiiSiSiSiifliii 
The combined effect of multi-path interference through 
several nearly degenerate levels of a single dot and corre- 
lations is discussed as one of the possible reasons for the 
puzzling behavior of the transmission phase measured in 
an Aharonov-Bohm interferometer li'* Without a direct 
hopping between the dots the setup in Fig. ^ E can ei- 
ther be viewed as two dots each having a single level or 
a single dot with two levels. For certain parameter sets 
the cases Fig. Q] D and E are equivalent (see below) 

The numerical renormalization group (NRG) is a re- 
liable method to investigate physical properties of sys- 
tems with local Coulomb correlations^^i^ and was very 
successfully used to study the linear transport through a 
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single-level dot^^'^i^ Unfortunately, the computational 
resources required strongly increase with the number of 
interacting degrees of freedom. This usually limits the 
applicability of NRG to mesoscopic systems with two 
spin-degenerate levels2Si22i21ii2ii2SiS and even in these 
cases a complete analysis of the parameter dependen- 
cies is practically impossible. For a linear chain of dots 
using special symmetries and a direct relation between 
the linear conductance and phase shifts it was possible 
to compute G at T = for up to four dotS:^^'^'^ 

Thus, alternative reliable methods are needed to 
systematically investigate parameter dependencies and 
study more complex systems with a larger number of 
degrees of freedom. Simple methods such as low-order 
perturbation theory or the self-consistent Hartree-Fock 
approximation (SCHFA) already fail to correctly describe 
the Kondo physics ^'^ in a single-level dot and are thus in- 
appropriate. Equations-of-motion techniques were suc- 
cessfully used to compute G{Vg) at high temperatures, 
but even sophisticated decoupling schemes only quali- 
tatively capture the small-T Kondo regimei32^ Similarly, 
perturbation theory in the coupling of the mesoscopic 
system to the leads works at large T, but fails for the 
large conductance resonances at small T.^^ In tempera- 
ture and parameter regimes in which charge fluctuations 
dominate a real-time renormalization group method^S. 
can be used;-^ A zero-temperature method that has fre- 
quently been used is exact diagonalization of a small clus- 
ter of lattice sites containing the dot followed by an em- 
bedding procedure!^ It has been criticized22i^ to pro- 
duce severe artifacts in G{Vg) and even in the case of 
a single-level dot, where at least the qualitative behav- 
ior comes out correctly, the quantitative agreement with 
NRG data is rather poor.**^ Also the mean-field slave- 
Boson approach'*'^ was used, first to study the single-level 
dot. Within this method G{Vg) falls off too quickly away 
from the plateau region as compared to the Bethe ansatz 
and NRG data^^*^ Furthermore, the method does not 
capture the correct dependence of G{Vg) on a magnetic 
field that lifts the spin degeneracy of the single leveli^ 
Despite these shortcomings the method was also used for 
more complex dot systems.-^ Recently, a method based 
on the Gunnarsson-Schonhammer variational approach^ 
has been applied to the dot with a spin-degenerate leve^ 



and a linear chain of three dots.^ It remains to be seen 
if this approach can with equal success be used in the 
presence of a magnetic field and for multi-path systems. 

We here propose a novel T — approach to 
study transport through mesoscopic systems with lo- 
cal Coulomb correlations that is based on the func- 
tional renormalization group (fRG). The fRG was re- 
cently introduced^ as a new powerful tool for study- 
ing interacting Fermi systems. It provides a systematic 
way of resumming competing instabilities^ and goes be- 
yond simple perturbation theory even in problems which 
are not plagued by infrared divergences iSiS The fRG 
procedure we use starts from an exact hierarchy of dif- 
ferential flow equations for the one-particle irreducible 
vertex functions^^Si^^iSiS i.e., the self-energy, the irre- 
ducible 2-particle interaction etc. It is derived by replac- 
ing the free propagator by a propagator depending on an 
infrared cutoff A and taking the derivative of the gener- 
ating functional with respect to A. We will apply two 
different truncation schemes to the models sketched in 
Fig. n Both result in a simple set of coupled differential 
equations that can easily be integrated numerically and 
in certain cases even analytically. By an explicit compar- 
ison of our results for G{Vg) and the occupancies to NRG 
and Bethe ansatz data we show that the method works 
well for surprisingly large U/T. The method can easily 
be extended to other setups and thus provides a useful 
and reliable tool to investigate transport through quan- 
tum dots. In two short publications it was earlier used 
to compute G{Vg) for a single dot with Luttinger-liquid 
leads^^^ and for a spin-polarized parallel double dotuS 

This paper is organized as follows. In Sec^Jwe present 
our general model. We then derive the fRG flow equa- 
tions in Sec. IIIII In Sees. II VI VIII we discuss our results 
for the gate voltage dependence of the linear conductance 
of the systems of Fig. ^ and for the cases shown in A-D 
we explicitly compare them to existing NRG and Bethe 
ansatz data. For each situation we judge the quality of 
our approximate method and discuss the limits of its ap- 
plicability. In Sees. IVIVIII fsetup of Fig. ^B, D, and E) 
we go beyond situations that have previously been stud- 
ied and present novel results. We conclude in Sec. IVIIII 
with a brief summary and an outlook. 
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II. THE MODEL 

Our general Hamiltonian 

H = -fflcad + ^^dircct + H^ot + H^on-p (1) 

consists of four parts. For simplicity the two leads are 
assumed to be equal and described by a tight-binding 
model 

oo 

i?icad = -r^^ ^ (c|„_^ ;C„+i,^,, + H.c.) (2) 

<y l=L,Rm=0 

with c'^^ I being the annihilation (creation) operator for 
an electron with spin direction a =t, i localized on lattice 
site TO of the left (/ = L) or right {I = R) lead. The 
hopping matrix element in the leads is denoted by r. To 
describe the system sketched in Fig.QlC we need a term 
that directly links the left and right lead 

-ffdiioct = ^ (tLR4),a,R^0,a,L + H.C.^ . (3) 

(7 

The coupled quantum dots are modeled as 

o" j>j' 

with annihilation (creation) operators dl'p^ and the site 

occupancy rij^^ = dj^dj^rj. With this dot Hamiltonian 
we equally describe multi-level and multi-dot situations 
or even mixtures of both. The dot or level index j runs 
from 1 to N. The level positions ej,a- consist of constant 
parts e° ^ and a variable gate voltage Vg that is used to 
shift the energy levels: ej^cr — ^'j „ + Vg- Our notation 
allows for both spin-degenerate levels and levels that are 
split by a magnetic field. For j's that belong to different 
dots it is meaningful to introduce nonzero hopping matrix 
elements tjji . The 2-particle interaction Uj'^ vanishes 
for j = j' and cr = a' . Besides this restriction we allow 
for both intra-level interactions of spin up and spin down 
electrons as well as inter-level interactions. Finally, the 
coupling between dot and lead states is given by 

^^coup = - E E + H.c.) (5) 

<^ l=L,R j 

with overlaps tj, that enter the energy scale of the dot 

level broadening F^ = 7r|tjppiead, where piead denotes 
the local density of states at the end of each semi-infinite 
lead. As usual we later take piead to be energy indepen- 
dent (wide-band limit). The choice of the nonzero t^- and 



tj,j> introduces a spatial geometry between the dots. For 
ring-like geometries (Fig. C and E) one might also be 
interested in studying the effect of a magnetic flux pierc- 
ing the ring. This can be modeled by including appro- 
priate phase factors in the hopping matrix elements tLR, 
t^-, tjj' along the ringiSSiSiiS All the systems shown in 
Fig. n follow from this Hamiltonian by special choices of 
the parameters. Obviously, it can also be used to model 
other cases of interest. 



III. THE METHOD 

Applying the Landauer-Biittiker formalism^ the T = 
linear conductance for noninteracting electrons can be 
expressed as 

GiVg)^G^iVg) + G^iVg) , (6) 

with 

GAVg)^j\UO,Vgf , (7) 
where the transmission is given by 

%{e, Vg) = 271x2 ga^^^L;o,aM£ + *0) (8) 

and the energy e is measured relative to the chemical po- 
tential /i. Here Q denotes the retarded one-particle Green 
function (the resolvent in the present single-particle prob- 
lem) of the system. For the leads described by Eq. Q 
the local density of states at site to = is (if they are 
decoupled from the rest of the system) 



Acad(£) = yiT^^TT+T^ ei2T-\e + fi\) . (9) 

The transmission phase a^iVg) is given by the phase of 

The noninteracting Green function matrix element in- 
volving lead states Go,a,L-o.cr,B. can be related to matrix el- 
ements involving only the dot states with quantum num- 
bers j,a and pioad- The explicit form of this relation 
depends on the choice of nonzero tj . It can be obtained 
using a standard projection technique^*^ and is based on 
the relation 

PGip{z)P = (10) 
zP - PHipP - PHipQ {zQ - QHipQ)-^ QHipP'] ' , 

where Hip is an arbitrary single-particle Hamiltonian, 
Gip{z) = [z — Hip] ^ the resolvent, and P, Q are projec- 
tors with P + Q = 1. 

For Fermi liquids, as discussed here, the same T = 
relations Eqs. I©-® between matrix elements of the one- 
particle Green function and the conductance hold in the 
presence of interactions!^ To obtain Go- and we thus 
have to determine the one-particle Green function of the 
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dot system in the presence of the interaction and the 
leads. From this also the level occupancies (nj,cr) can be 
determined by integrating Gj,a-.j,cr{z) along the imaginary 



axis.--^ 

As a first step we integrate out the noninteracting leads 
within a functional integral representation of our many- 
body problem.^- They provide a frequency dependent 
one-particle potential for those indices j for which 7^ 0. 
On the imaginary frequency axis it is given by 



I 



(11) 



where gi^g^^^icu) denotes the spin- independent Green func- 
tion of the isolated semi-infinite leads taken at the last 
lattice site 



Dlcad 



2r2 



1- Wl- 



4t2 



{iuj + fj,y 



(12) 



After this step, instead of dealing with an infinite system 
we only have to consider the dot system of interacting 
levels. We here take /i = 0. In the computation of the 
Green function projected onto the dot system the sum of 
the dot Hamiltonian at Uj'^ — and Vj°p'^{ioj) can be 
interpreted as a frequency dependent "effective Hamilto- 
nian" and in the following will be denoted by h{iuj). 

We now set up our fRG scheme. As a first step we 
replace the noninteracting dot propagator Qq, as obtained 
from Eq. (0J and the projection of the leads, by a cutofl"- 
dependent propagator that suppresses the infrared modes 
with Matsubara frequency |a;| < A, 

g^{iLu) = e(|c^| - A)5o(«^) = e(|w| - a) [iuj - h{iLu)]-^ 

(13) 

with A running from 00 down to 0. For convenience we 
use a sharp cutoff (see below). Using Qq in the gener- 
ating functional of the irreducible vertex functionsS^ and 
taking the derivative with respect to A one can derive an 
exact, infinite hierarchy of coupled differential equations 
for vertex functions, such as the self-energy and the irre- 
ducible 2-particle interaction. In particular, the flow of 
the self-energy (1-particle vertex) is determined by 
T,^ and the 2-particle vertex F^, while the flow of F^ is 
determined by S^, F^, and the flowing 3-particle vertex 
Fg . The latter could be computed from a flow equation 
involving the 4-particle vertex, and so on. At the end 
of the fRG flow E^^^^ is the self-energy S of the orig- 
inal, cutoff-free problem we are interested imSi^SiS^ A 
detailed derivation of the fRG flow equations for a general 
quantum many-body problem that only requires a basic 
knowledge of the functional integral approach to many- 
particle physics^ and the application of the method for 
a simple toy problem are presented in Ref. J3. In practi- 
cal applications the hierarchy of flow equations has to be 
truncated and S^^*^ only provides an approximation for 
the exact S. As a first approximation we here neglect the 
3-particle vertex (irreducible 3-particle interaction). The 
contribution of Fg to F'^ is small as long as F'^ is small. 



because F^ is initially (at A = 00) zero and is generated 
only from terms of third order in F^. Furthermore, F^ 

stays small for all A if its initial values C/J^T (for the 
details of the initial conditions, see below) are not too 
large. By explicit comparison to NRG data below we 
will clarify the meaning of "not-too-large" in the cases 
of interest. This approximation leads to a closed set of 
equations for F'^ and given by 



d_ 
dX 



1 

2^ ^ 

2,2' 



5^(2,2') F^(l',2';l,2) 



and 



(14) 



Af^I', 2'; 1, 2) = ^ 5] ^ (?^3, 3') S\A, 4') 

3,3' 4,4' 

X F^(1',2';3,4)F^(3',4';1,2) 

- F^(l', 4'; 1, 3) F^(3', 2'; 4, 2) - (3 ^ 4, 3' ^ 4') 

+ F^(2', 4'; 1, 3) F'^(3', l'; 4, 2) + (3 ^ 4, 3' ^ 4' 



.(15) 



The labels 1, 2, etc., are a shorthand for the quantum 
numbers of the one-particle basis j, a and the Matsubara 
frequencies and the summation stands for a sum over the 
quantum numbers and an integral over the frequency. 
The full propagator is given by the Dyson equation 



and the so-called single-scale propagator by 



(16) 



(17) 



The order of the projection onto the system of dots and 
the introduction of a cutoff can equally be interchanged, 
leading to the same set of fiow equations. 

We next implement our second approximation: the 
frequency-dependent flow of the renormalized 2-particle 
vertex F^ is replaced by its value at vanishing (exter- 
nal) frequencies, such that F^ remains frequency inde- 
pendent. Since the bare interaction is frequency indepen- 
dent, neglecting the frequency dependence leads to errors 
only at second order (in the interaction strength) for the 
self-energy, and at third order for the vertex function at 
zero frequency. As a consequence, also the self-energy 
becomes frequency independent. Then ej,a- + ^j^a-;j,a- can 
be viewed as the effective (single-particle) level position 
and tj,j' + ^j^a;j',(7, with j ^ j' as the effective (single- 
particle) hopping between the levels. In general both de- 
pend on all parameters of the problem, in particular the 
two-particle interaction. This interpretation will later be 
helpful to gain further insights in our results. 

For a sharp frequency cutoff the frequency integrals on 
the right-hand side (rhs) of the flow equations (|14|l and 
(|15|l can be carried out analytically. At this point one has 
to deal with products of delta functions S{\uj\ — A) and 



5 



expressions involving step functions Q{\uj\ — A). These 
at first sight ambiguous expressions are well defined and 
unique if the sharp cutoff is implemented as a limit of in- 
creasingly sharp broadened cutoff functions O^, with the 
broadening parameter e tending to zero. The expressions 
can then be conveniently evaluated by using the follow- 
ing relation, valid for arbitrary continuous functions /: 

S,{x - A) f[QAx - A)] ^ S{x - A) f f{t)dt, (18) 

Jq 

I 



where 5e = Q'^- Note that the functional form of for 
finite e does not affect the result in the limit e ^ 0. For 
the approximate flow equations we then obtain 

ui = ±A 2,2' 

and 



|^r^',2';l,2 = ^ E E E GtA'H^) r^',2';3.4r^',4';1.2 + ^3^ (*^) ^^4' M 

uj=±A 3,3' 4,4' 

^ [^ri'.4';l,3 r3',2';4,2 + r2',4';l,3 r3',l';4,2] | (20) 



where the lower indexes 1, 2, etc. now stand for the single- 
particle quantum numbers j, a (not frequencies) and 

Q^{ioj)^{Q^\iuj)-Y.^Y^ . (21) 

Note that in contrast to C/^ in Eq. (|16|l . Q'^ contains the 
cutoff-free noninteracting propagator. Thus, by taking a 
sharp cutoff the explicit cutoff function completely disap- 
pears from the flow equation, making it smooth and easy 
to integrate numericallyiS^ At the initial cutoff A = oo 
the flowing 2-particle vertex Vy 2'-i 2 given by the an- 
tisymmetrized interaction (see below). Equation IjlSfl 
preserves the antisymmetry under the exchange of the 
first and second two indices. This reduces the number of 
independent F^, j/.i 2- Furthermore, Eq. (|20|l preserves 
the spin conservation of the Hamiltonian. One can ei- 
ther use these symmetries to reduce the number of both, 
flow equations and independent summations on the rhs 
of Eq. H2U|) , or directly implement this equation for sim- 
plicity. In complex cases the first will considerably speed 
up the numerical integration of the flow equations. To 
exemplify this we consider the single-level case. Then 
the only physical interaction is the one between spin up 
and spin down electrons and Eq. (^01 reduces to a single 
equation with no sum over quantum numbers left on the 
rhs. 

The flow is determined uniquely by the differential flow 
equations and the initial conditions at A = oo. The 
flow of the 2-particle vertex starts from the antisym- 
metrized bare 2-particle interaction?^ -^i',2';i.2 while m- 
particle vertices of higher order vanish in the absence of 
bare m-body interactions with m > 2. The self-energy 
at A = oo is given by the bare one-particle potential, 
that is, by those one-particle terms which are not in- 
cluded already in Qq . For fixed quantum numbers j, a in 

Eq. Q these are terms of the form — J^a' j' ^j 'f "'^^ 
a numerical solution the flow starts at some large finite 



initial cutoff Aq. Here one has to take into account that, 
due to the slow decay of the rhs of the flow equation for 
at large A, the integration of the flow from A = oo 
to A = Aq yields a contribution which does not vanish 
in the limit Aq — s- oo, but rather tends to a finite con- 
stant. Since ^^2'(*^) ~* <52.2'/(*w) for |cj| = A — s- oo, this 
constant is easily determined as 

-J- hm TdA 5^e-°^^/,,2';,2 

= ^E^i'^2^i.2. (22) 

2 

This exactly cancels the bare one-particle term (see 
above). In summary, the initial conditions for the self- 
energy and the 2-particle vertex at A = Ag ^ oo are 

T,^"y = (23) 

l'l'°2';l,2 = ^l',2';l,2 ■ (24) 

For the flow at A < Aq the factor e*'^°^ in Eq. (^HJ can 
be discarded. 

In the following sections we use the above fRG-based 
approximation scheme to study the conductance, the 
transmission phase, and the level occupancies for the 
geometries sketched in Fig. ^ For the single-dot case 
we also apply an even simpler scheme for comparison, 
where the flowing 2-particle vertex is replaced by its ini- 
tial value. The system of flow equations then reduces to 
Eq. lEl). 

IV. SINGLE DOTS 

As a first example we study the single dot case of Fig.Q] 
A. Since = 1 we suppress the dot index j = 1. Fur- 
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thermore, the Hamiltonian conserves the spin and the 
self-energy and Green function is diagonal in the spin in- 
dex a. The same holds for the other setups studied. The 
level energies are e-f = Vg — ?i/2, ei = Vg+H/2, where H 
denotes a magnetic field, and the interaction between up 
and down electrons occupying the dot level is U. We here 
consider the case with left-right symmetry of the level- 
lead hoppings = — t' , but the results can easily 
be generalized. Note that due to the shift of the density 
by —1/2 in the last hne of Eq. = corresponds 

to half filling of the dot. The projected noninteracting 
propagator is given bjii^ 



1 



iuj - {Vg + (jH/2) + irsgn(cj) ' 



(25) 



where on the rhs a = t = +1 and a = i = —1. We have 
performed the wide-band limit that leads to an energy- 
independent hybridization F = F'^ + F''*^ = 2'Kt'^picad- 
This is achieved by replacing t rjT and t' y/rjt' and 
taking 77 — s- cx). The flow equations for the effective level 
positions = Vg+an/2 + arc 



d_ 
dA 



2tt 



(A + F)2 + (F/)2 



= V„ 



(26) 



(tH/2 and a 



with the initial condition \M=°° 
denoting the complement of a. The cutoff-dependent 
propagator Q^{iu!) follows from Qo^cr{ii^) by replacing 
Vg + aH/2 V^. As already mentioned above, us- 
ing symmetries the flow of the 2-particle vertex can be 
reduced to a single equation for the effective interaction 
between spin up and spin down electrons U^. It is de- 
fined by = T^g.^g, has the flow equation 



d_ 



^ [^f (*^) GtH^) + Gt(^^) Gii^^) 



2{U^)'v{^V,^/n 



(A + F)2 + (1M)2 (A + F)2 -I- (FA)2 



(27) 



and initial condition U^=°° = JJ . 

Within this approximation the dot spectral function at 
the end of the fRG flow is given by 



1 



p^(w) = - 



F 



^ (c^-V;)^ + F2 



(28) 



with K = V^=^, that is a Lorentzian of full width 2F 
and height l/(7rF) centered around V„. For the single 
dot, instead of using Eq. (jHJ the transmission can be ex- 
pressed in terms of the dot spectral weight at the chem- 
ical potential^ 



Ga{Vg 



(29) 




FIG. 2: (Color online) Gate Voltage dependence of the ef- 
fective level position and interaction for a single-level dot at 
U/V = A-K and H = 0. Dashed line: V = V^=° without flow 
of the interaction. Dashed-dotted line: V = V^^'^ with flow 
of the interaction. Solid hne: effective interaction ?7"^^°. 
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FIG. 3: (Color online) Gate Voltage dependence of the linear 
conductance for the same parameters as in Fig.|21 Solid line: 
exact Bethe ansatz solution from Ref. I21I . Dashed line: fRG 
approximation without flow of the interaction. Dashed-dotted 
hne: fRG approximation with flow of vertex. 



Within our approximation the spectral weight and thus 
the conductance directly follows from the effective level 
position Va at the end of the fRG flow. 

In the following we consider strong couplings C//F 3> 1 
and start out with the case H = 0. The latter implies 
that = Vf" = and Gt = G^. The spectral func- 
tion computed at Vg = using NRG (that is beheved to 
be very close to the exact spectral function) shows a sharp 
Kondo resonance of height l/(7rF) and width Tk, where 
Tk is the Kondo temperature, located around uj = 0. It 
has two additional broader features at w = ztU/2 (Hub- 
bard bands). For -U/2 ^ Vg ^ U/2 the Kondo peak is 
pinned at (close to) the chemical potential and has fixed 
height. ^'^ Using Eq. this leads to the plateau of width 
U in G{Vg) discussed in the introduction. Although our 
approximate spectral function neither shows the narrow 
Kondo resonance nor the Hubbard bands, we show that 
by a pinning of Va- it captures the pinning of the spectral 
weight. We reproduce the line shape of the conductance 
quantitatively up to very large U/T. 

To gain some analytical insight we first replace — > 
U and solve Eq. (|26|) . The solution is obtained in implicit 
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form 



vJijv) - jJojv) _ Jojvg) 
vY,{v) - jYoiv) ~ Yoivg) 



(30) 



with V = Vtt/U, Vg = VgTr/U, 7 = Tn/U, and Bessel 
functions J„, F„. For \Vg\ < Vc, with Vc = Vcn/U being 
the first zero of Jq, i.e., Vc = 0.7655 U, this equation has 
a solution with a smaU \V\. For [/ 3> F the crossover to 
a sohition with \V\ being of order U (for \Vg\ > Vc) is 
fairly sharp. Expanding both sides of Eq. for small 
\v\ and \vg\ gives V = Vg exp [— [//(ttF)]. Inserting this in 
Eq. (|28|l leads to an exponential pinning of the spectral 
weight at fj.. In Fig. |2 we show V{Vg) for U/T = An 
as the dashed line. For |Fg| » U/2 we find V ^ Vg - 
sgn(V^) U/2. Applying Eq. the exponential pinning 
of the spectral weight at /i = for small \Vg\ and the 
sharp crossover to a y of order U when \Vg\ > Vc leads to 
the plateau in G{Vg) shown as the dashed line in Fig. O 
For [/ > r the width of the plateau is 2Vc = 1.531C/, 
which is larger than the width U found from the exact 
Bethe ansatz solution'^ ^ (solid line in Fig. O and with 
NRG. 

This can systematically be improved including the flow 
of the interaction Eq. (|27|l . As the dashed-dotted and 
solid lines in Fig. |2] we show the effective level position 
and the effective interaction at the end of the fRG flow. 
The dashed-dotted line in Fig. Ois the resulting G{Vg). 
Away from Vg — the effective interaction first decreases 
and the plateau becomes narrower. The approximate 
conductance then agrees quantitatively with the exact 
Bethe ansatz result (solid line in Fig.O. For decreasing 
U/r the agreement of fRG and Bethe ansatz data sys- 
tematically improves and for U/T 6 the two curves can 
barely be distinguished. The agreement becomes only 
slightly worse than in Fig. Ofor U/T as large as 25, the 
largest U /F for which Bethe ansatz data are available in 
the literatureiSi Also for more complex dot geometries 
including the flow of the effective interaction consider- 
ably improves the agreement with NRG and from now 
on we exclusively use this truncation scheme. 

We next consider the case of finite magnetic fields. For 
H > the Kondo resonance in the NRG solution of the 
spectral function splits in two peaks with a dip ai lu — 0, 
resulting in a dip of G{Vg) at Vg ^ 0. In Figs. ^ and 
13 we compare the total G = Gf -I- Gj^ and partial G| 
conductance obtained from our fRG truncation scheme 
including the flow of the interaction and from NRG^^ for 
U/T = Stt and different Ti.. In the caption we give Ti. in 
units of r^i?G _ o.llGF, where Tj^^^ is determined from 
the width of the Kondo resonance at the particle-hole 
symmetric point = using NRGi^ The agreement 
between NRG and fRG results is excellent. In particular, 
at Vg = the NRG and fRG data for Ti. = Tj^^^ are 
almost indistinguishable. 

This observation suggest that our approximate scheme 
indeed contains the Kondo temperature Tr- that de- 
pends exponentially on a combination of U and the level 




FIG. 4: (Color online) Gate Voltage dependence of the total 
conductance G of a single dot with U/T = 3n and Ti/T = 0, 
0.058, 0.116, 0.58 from top to bottom. In units of the Vg ^ 
Kondo temperature TJ^^'^/F = 0.116 these fields correspond 
to H = 0, 0.5Tk^°, Tk^°, and ST]}''^ . Solid line: NRG data 
from Ref. |6^ . Dashed line: fRG approximation with flow of 
vertex. 




FIG. 5: (Color online) Gate Voltage dependence of the partial 
conductance G| of a single dot for the same parameters as in 
Fig.Bl 




FIG. 6: (Color online) The Kondo scale Tk as a function of 
U for different Vg-. Vg = (solid), Vg/T = 1 (short dashed), 
Vg/r = 2 (dashed-dotted), and Vg/T = 4 (long dashed). In- 
set: The fitting parameter c as a function of (Vg/T)'^. 
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(31) 



The prefactor of the exponential depends on the details 
of the model considered. To leading order its U and Vg 
dependence can be neglected. To investigate the appear- 
ance of an exponentially small energy scale within our 
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approximation we define a Kondo scale Tk{U, Vg) by the 
magnetic field required to suppress the total conductance 



down to half the unitary limit G{U, Vg,Ti. — Tk) 



For fixed U '3> T this definition is meaningful for gate 
voltages which for Ti. — are in the conductance plateau. 
In Fig.Elwe show Tk{U, Vg) for different Vg as a function 
of J7 on a linear-log scale. The curves can be fitted by a 
function of the form [see Eq. l|lJT|l ] 



/(C//r) = aexp 



U 



'U 



with Vg dependent coefficients a, b, and c. On the scale 
of Fig. 1^1 the original data and the fits are indistinguish- 
able. We find that b{Vg) barely changes with Vg and 
is given by 6 « 0.32, which is close to the exact value 
7r/8 « 0.39 [see Eq. (|?T]| ]. Furthermore, the prefactor a 
depends only weakly on Vg and c(Vg) increases approxi- 
mately quadratically with Vg /T as shown in the inset of 
Fig. El Both these results are consistent with the behav- 
ior of the exact Kondo temperature Eq. (|31|l . We thus 
conclude that TxiU^Vg) can be estimated from the Ti 
dependence of G obtained within the fRG. The Kondo 
temperature can also be obtained from the local spin 
susceptibilityii^ Computing this using our fRG scheme 
leads to results equivalent to the ones discussed above. 

For the single dot at T = the exact conductance, 
transmission phase, and dot occupancies are directly re- 
lated by a generalized Friedel sum rulei? Ga/{e^/h) = 
sin^ (tt (ricr)), a„ — n {n„) . As < (71^) < 1 the ar- 
gument of sin^ is restricted to a single period and the 
relation between G^, {n,y), and is unique. In many 
approximation schemes the Friedel sum rule does not 
hold exactly. Within our method we map the many-body 
problem onto an effective single-particle one for which the 
Friedel sum rule is fulfilled. For gate voltages within the 
Ti = conductance plateau the (spin independent) dot 
fiUing is 1/2 and the (spin independent) phase is 7r/2. 
For sufficiently large U /F the crossover to (no-) = 1 and 
tto- = TT to the left of the plateau as well as {ua-) = and 
tto- = to the right is fairly sharp. 

As a more complex single-dot problem we study the 
case of a dot embedded between two leads, that are also 
coupled directly (see Fig. ^ C) . In that case the conduc- 
tance shows the characteristics of both the Kondo and 
the Fano effect, as was discussed earlier based on NRG 
results J£ We here focus onH — 0. 

To derive expressions for the dot Green function and 
the conductance we first consider a three-site system of 
the dot and the last sites of the left and right leads. The 
effect of the other lead sites is taken into account by 
projecting them out, which leads to a single-particle po- 
tential on the last sites of the leads similar to Eq. Hll() 
but with t^j replaced by r. At /7 = the projected "effec- 
tive Hamiltonian" [see Eq. H10|l ] for both spin directions 
reads (for the hopping between the leads and dot we still 




FIG. 7: (Color online) Gate Voltage dependence of the total 
conductance G of a single dot in the Fano geometry of Fig. Q 
C with U/r ^ 5.057r, H ^ 0, and different Tlr. Symbols: 
NRG data from Ref. i2Q (circles Tlr — 0, squares Tlr — 0.1, 
diamonds Tlr = 0.6, triangles Tlr — 1). Solid lines: fRG 
approximation with flow of vertex. 



assume = 



T^glcad(«^) 
tLR 
-t' 



tLR 



-t' 



-t' 



(32) 



Va 



in a basis left lead — right lead — dot. We here consider 
t'jtLR G R (for the model with an enclosed magnetic flux, 
see below). The Green function matrix element entering 
Eq. ||HJ) is given by the 1 — 2 matrix element of the resol- 
vent of ha{iuj), Eq. (|^ . Within our truncation scheme 
the same holds for [/ > if is replaced by the effec- 
tive dot level position V^^^ in /i3^3(ia;). To determine 
the flow equations for V^ and we employ Eq. t\10\\ 
a second time and project onto the dot site. After tak- 
ing the wide-band limit the cutoff-dependent dot Green 
function is 



V' 



tLR + tsgn(cj) 



(33) 



with tifl — Trtijjpioad- The first lines of Eqs. (|26|) and 
()27|l also hold in the presence of the direct hopping t^n, 
which completes the derivation of the flow equations. 

At J7 = the interference of the direct energy indepen- 
dent transmission channel and the resonant channel leads 
to a Fano line shape of G{Vg)^ The line shape result- 
ing from the presence of both the Kondo and the Fano 
effect is shown in Fig.[7|for^& U/T « 5.057r and different 
strength of the direct hopping between the left and right 
lead. The latter is measured by the U = transmission 
probability TLR{tLR) that would result if only the direct 
link would be present, that is for t' = 0. For Tlr = 
one recovers the single-dot problem discussed above. In 
the opposite limit of perfect direct transmission Tlr. = 1 
the U — Fano anti-resonance at = (see the g = 
curve in Fig. 1 of Ref. I25I) is extended to a broad region 
with G{Vg) « due to the presence of the Kondo effect. 
The results for < Tj^r < 1 lie in between these two 
limiting cases. The fRG data again agree well with the 
NRG resuhs. 
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Similar to the single dot without a direct link for the 
present geometry one can derive unique relations be- 
tween the conductance, the transmission phase, and the 
occupancyiiS& Thus the last two observables do not carry 
any additional information and we do not show them 
here. In addition to the gate voltage dependence also 
the dependence on an enclosed magnetic flux was stud- 
ied. Replacing ^ exp t^ji a flux can easily be 
included in our approach and we reproduce the results of 
Ref.ii 

Our results for the single-level dot show that although 
the fRG truncation scheme is set up for small U, it works 
quantitatively for very large U. For 7i = the method 
captures the aspect of Kondo physics essential to ob- 
tain the plateau in GiVg), namely the pinning of spectral 
weight. From this one can expect that the method gives 
reliable results also for more complex dot systems, as will 
be explored below. Using a truncation scheme in which 
the full frequency dependence of the 2-particle vertex is 
kept (leading to a frequency dependent self-energy) it was 
shown that one can also reproduce the Kondo resonance 
and Hubbard bands of the spectral function^^ although 
with a much higher numerical effort S 

V. LINEAR CHAINS OF DOTS 

We next study a linear chain of N single-level dots as 
shown in Fig. ^B. For this problem one already has to 
specify a variety of parameters. The inter-dot hopping 
matrix elements might not only be restricted to nearest 
neighbors but extend over a longer range. Depending 
on the experimental realization of the chain in addition 
the hoppings of equal range might not all be the same 
but explicitly depend on the pair of dots considered. Be- 
sides the on-site interaction (that might depend on j) 
one could introduce longer-range interactions. The on- 
site energies Sj^a- might vary from site to site. Within our 
approach all these situations can be investigated, as the 
parameters only enter the projected noninteracting prop- 
agator and the initial conditions for the 2-particle vertex. 
We here consider the cases = 3 and = 4, but also 
A'' = 10 has been studied with our method. We focus 
on 7i = 0, nearest-neighbor inter-dot hoppings of equal 
amplitude, equal on-site energies, and equal local inter- 
actions. We verified that including a nearest-neighbor 
interaction and a small variation of the inter-dot hop- 
pings only weakly affects our results. 

At U = the projected "effective Hamiltonian" is a 
tridiagonal N x N matrix with off-diagonal entries 

ha--j,j±iiiuj) = ~t (34) 

and diagonal elements 

h^,jj (luj) ^Vg-i sgn(cj) (rf 5,- 1 + T^Sj^n) ■ (35) 

The (spin-independent) projected noninteracting propa- 
gator then follows from the matrix inversion 

Go,a{ii^) ^ [ii^ - ha{iu!)r^ ■ (36) 




FIG. 8: (Color online) Gate Voltage dependence of the con- 
ductance G and the total occupancy (n) of a linear chain of 
= 4 spin-degenerate dots (see Fig.^BVwith fZ/F — 2-k and 
t/r = 1. Circles: NRG data from Refj^l Solid line: £RG 
approximation with flow of vertex. 

The flow equations for the present setup are obtained 
by inserting this into Eq. H21|l and using the result- 
ing cutoff-dependent full propagator in Eqs. ifT^ll and 
(I20|l . The initial conditions are given by Eqs. 123|l and 
(HHl. Although we start out with a purely local inter- 
action, during the fRG flow Eq. H2U|) all types of two- 
particle interactions not forbidden by spin conservation 
and antisymmetry of the vertex, e.g., long-range density- 
density interactions and pair hoppings, are generated in 
F^, , , / . We observe the same for the setups 

of Figs. 1 D and E. With this 2-particle vertex the flowing 
self-energy Eq. H19|) becomes a full NxN matrix. During 
the fRG flow the interaction does not only lead to a renor- 
malization of the level positions and the nearest-neighbor 
hoppings, but also generates longer-range single-particle 
hoppings. Applying Eqs. l|H|)- (|l()|l the conductance can be 
computed from the 1 — matrix element of the projected 
Green function at A = by (note that = 0^=°) 

GAVg) = ^4TfT^\g^=,%{0)\' . (37) 

The occupancy of dot j follows from integrating t/j^~°(iw) 
over Lu. 

For the special case of equal nearest-neighbor inter- 
dot hoppings t, equal hybridizations Ff ~ Trjif'ppicad = 
F^ = 7r|f^ppioad = r/2, equal local interactions U, and 
equal Sj^^ = Vg, G{Vg) at T = was computed using 
NRGiSSiS^ For this class of parameters the Hamiltonian 
has a high symmetry and the conductance is related to 
scattering phase shifts. The latter can be extracted from 
energies of the system which allowed to obtain results for 
a chain of up to four dots. In Fig. |Slwe show a compari- 
son of the NRG and fRG results (for G and the total dot 
occupancy (n) = ^ (nj>)) for A^ = 4, U/T = 27r, and 
t/F = 1. For gate voltages at which the four-site chain 
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is occupied by an odd number of electrons one finds a 
transmission peak of unitary height. Each corresponds 
to transport through mainly one of the four dot levels. 
For U /T — 2ti the plateau-like line shape of the indi- 
vidual resonances induced by the Kondo effect is only 
weakly developed. It becomes more pronounced if larger 
U /T are considered. The overlap of the resonances de- 
creases for increasing t/T. For even TV the width of all 
the resonances is almost equal and of order U /N . In 
a wide region around = 0, GiVg) is small because 
close to half-filling of the chain charge fluctuations are 
strongly suppressed. For odd N , as in Fig.^ half-filling 
corresponds to an odd number of electrons on the chain. 
This implies a wide plateau of unitary height around Vg 
due to the suppressed charge fluctuations (compare the 
solid lines in the upper U /T = tt and lower panel U — 
of Fig. ISJ and only the resonances away from Vg — Q 
have width U /Nw^ The (spin-independent) transmission 
phase ac{Vg) and the total occupancy are related by 
= 2acr/7r>2iS Across every transmission resonance 
the phase continuously changes by tt in analogy to the 
behavior of the phase at J7 = 0. Note that in contrast 
to the single-dot case a„ is not restricted to the interval 
[0,^]. 

For U /T smaller than the one of Fig. |Slthe agreement 
between NRG and fRG further increases, e.g. for iV = 4, 
U /T = 1.677r, and t/T — 4.167 both data sets are al- 
most indistinguishable (for the NRG data see Fig. 5 of 
Ref. |2^ . We encounter serious deviations of the fRG 
data and NRG results if U /V is significantly increased 
(i.e. for U/T = 8.33; for NRG data see Fig. 4 of Ref.lH). 
In that case regions of gate voltages appear in which some 
of the components of the effective interaction become or- 
ders of magnitude larger than the initial U leading to 
a conductance that is too small. For these Vg it is no 
longer justified to neglect the frequency dependence of 
the 2-particle vertex and the higher-order vertices. In 
particular, this happens for Vg on-resonance and in addi- 
tion for even N in the deep valley around Vg — Q. Indi- 
cations of this breakdown of our approximation scheme 
for U/T ^ 10 can already be observed in the Vg = Q 
valley, and the large | Vg \ resonances in Fig. in these 
gate- voltage regimes our approximation significantly un- 
derestimates G. 

In Fig. 1^1 we show how an increasing left-right asym- 
metry of the rf and affects the conductance for 
iV = 3 for small to intermediate t/T. The parameters 
are U/T = vr, with fixed F = Ff + Ff , and t/F = 1. For 
comparison in the lower panel we show the U — Q case. 
With increasing asymmetry the height of the peaks de- 
crease. While for fixed asymmetry the interaction leads 
to an enhancement of the outer two peaks, it suppresses 
the peak located around Vg = 0. Therefore, at J7 > 
the central peak disappears quickly if Ff/F^ increases. 
It is important to note that with increasing asymmetry 
the gate voltage dependence of the total occupation (n) 
and the transmission phase a„ (which are still related 
by {n) = 2aa/T^) barely change. Thus these observables 
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FIG. 9: (Color online) Gate Voltage dependence of the con- 
ductance G of a linear chain of A'^ = 3 spin-degenerate dots 
with asymmetric coupling to the leads. The parameters are 
U/T = TT (upper panel), U = Q (lower panel), t/V — 1, 
rf/rf = 1 (solid line), Tf /T^ = 2 (short dashed hne), 
Fi/rf = 4 (dashed-dotted line), and T^ /Vi = 9 (long 
dashed line). Note that here Vg is shown in units of F. 



do not show any indication of the disappearance of the 
central resonance. The same holds for the (effective) dot 
level-positions of the chain obtained by diagonalizing the 
"effective Hamiltonian" at the end of the fRG flow (com- 
puted in the presence of the leads) after disconnecting the 
three sites from the leads. In particular, for all Ff/Ff- 
studied one of the level energies (eigenenergies of the iso- 
lated chain at A = 0) goes through zero at = 0. For 
levels which are not too strongly (sufficient separation of 
the resonances) and not too asymmetrically coupled to 
the leads one would thus expect a peak at Vg = (see 
e.g. the U — case). Surprisingly, during the fRG flow 
the effective asymmetry of the coupling to the outer two 
levels (as obtained from the eigenvectors of the isolated 
chain at A = 0) is reduced compared to the initial value, 
while the asymmetry of the central level is significantly 
increased by the presence of the interaction. Further- 
more, the sum of the effective couplings of the central 
level to the left and right leads is increased which im- 
plies that the level is broadened. The combined effect of 
the asymmetry and broadening explains the vanishing of 
the central peak. For increasing peak separation, that is 
for increasing t/T, the central resonance becomes more 
robust against a left-right asymmetry. 

The above observation might be important for the in- 
terpretation of future experiments on linear chains of 
quantum dots because asymmetric dot-lead couplings are 
generic in realistic dot systems. As shown, such an asym- 
metry might lead to a two peak structure although a 
chain of three sites is studied. 
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VI. SIDE-COUPLED DOUBLE DOTS 

The side-coupled double dot of Fig. ^ D shows a very 
interesting low-energy physics. Close to half-filling it is 
either dominated by a two-stage Kondo effect for small 
tl 2 or the formation of a local "molecular" spin-singlet 
between the electrons on the two dots for large ti 2- Both 
regimes can be distinguished by the temperature depen- 
dence of the G(Va ) curves as is nicely discussed in Ref.l28l 
(see also Ref. 29) based on NRG data. In the two-stage 
Kondo effect, first at high energies the spin of the em- 
bedded dot gets screened at the Kondo scale of the 
system with ti^2 = 0. In the second stage the heavy 
quasi-particles screen the spin on the side-coupled dot at 
a much lower scale Tq. 

The U — projected "effective Hamiltonian" for spin 
direction cr reads (in the wide-band limit) 



+ rf)sgn(c^) -ii,2 

tl,2 £2,0- 



(38) 



As for the linear chain of dots the projected noninteract- 
ing propagator QQ{iuj) entering the flow equations H19() 
and (|20|1 follows by the matrix inversion Eq. H36|l . Start- 
ing out with a local interaction U (and maybe a nearest- 
neighbor density-density interaction U') all two-particle 
interactions not excluded by symmetries are generated 
(see above). Also in the self-energy matrix all four ma- 
trix elements flow, that is, both the level positions and 
the inter-dot hopping arc rcnormalized. As for the sin- 
gle dot we can use Eq. to compute the conductance 
from the spectral weight of the embedded dot with j = 1 
taken at A = Oi^ 

In Figs.lTUlandlTTIwe compare fRG results for G{Vg) to 
NRG data. Two different parameter sets with purely lo- 
cal interaction J7, = Ff = F/2, and ei.cr = £2,0- = Vg 
are used. In Fig. ^| the parameters are in the spin- 
singlet regime: U /T = 8, ti,2/F = AM The fRG and 
NRG data agree well. From the gate voltage dependence 
of the many-body eigenvalues of the double dot system 
without the leada^ it becomes clear that each of the 
plateau-like resonances centered around ±VJ can be un- 
derstood as the resonance in the single dot case, with the 
difference that the approximate width of the two peaks 



is A = U/2 + 2t 



1,2 



(C//4ti,2)2 



instead of U . 



Around Vg — both dots are half-filled and a local spin- 
singlet is formed leading to a suppression of the conduc- 
tance. The parameters U /T = 2 and ti,2/r = 0.2 of 
Fig. ^2 are taken from the two-stage Kondo regime, 
As U /F is moderate the first stage of the Kondo effect is 
not well pronounced which explains that the conductance 
away from half-filling (that is for \Vg\ > U/2) vanishes 
only slowly with increasing \ Vg\. For these gate voltages 
the fRG and NRG data are almost indistinguishable. The 
second stage of the Kondo effect generates a dip at a; = 
in the spectral function of dot j = 1 at half filling,^* that 
leads to the suppression of G near Vg = 0. For small 
ti^2 the scale C//ti 2 relevant in the second stage of the 




FIG. 10: (Color online) Gate Voltage dependence of the 
total conductance G and total occupancies (nj) of a spin- 
degenerate side-coupled double-dot system (see Fig.^D) with 
U/r = 8 and ^1,2/!" = 4. Solid line with left y-axis: NRG data 
for G from Ref. l28L Dashed line with left y-axis: fRG approx- 
imation of G with flow of vertex. Thick solid line with right 
y-axis: fRG approximation of (ni) (embedded dot). Thin 
solid line with right y-axis: fRG approximation of {n2) (at- 
tached dot). 




FIG. 11: (Color online) As in Fig. [Ullbut for U/F = 2 and 
ti,2/r = 0.2. 
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FIG. 12: (Color online) Gate Voltage dependence of the (spin- 
independent) transmission phase Ocr for the two parameter 
sets of Figs, irni and im 



Kondo effect becomes very large and, not surprisingly, 
our approximation deviates quantitatively (not qualita- 
tively) from the NRG data for |Vg| < U/2. For these 
gate voltages some of the flowing interactions become 
large and our truncation scheme is questionable. Within 
our approximation the parameter regime of very large 
U/ti 2, in which in addition to the two-stage Kondo ef- 
fect a Fano-like feature in G{Vg) was observedjiSS is out 
of reach. 

In Figs. llUl and llll in addition to G we show the dot oc- 
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cupancies {uj), with Uj = rij^-^ + rij^i, computed by fRG. 
Because of the large inter-dot hopping in Fig.^|the occu- 
pancies behave similarly for all V^. For small ti.2 (Tie. 111!) 
the side-coupled dot j = 2 is coupled very weakly to the 
leads and the particle number of this dot changes very 
abruptly across the transmission resonances. Close to 
Vg = the occupancies (rij) depend non-monotonically 
on Vg. This follows from an interaction-induced non- 
monotonicity of the effective level positions. Similar be- 
havior was observed in a spin-polarized parallel double- 
dot modeliSiS2i2£ In contrast the total dot charge (n) — 
(rii) -|- (712) is monotonic. Similar to the single-dot case 
it is directly related to the conductances^ 



G/{e^/h) = 2sin2((n) 7r/2) . 



(39) 



The (spin-independent) transmission phase for the pa- 
rameter sets of Figs. llOl andllllis shown in Fig.^| Exper- 
imentally this phase is accessible if the side-coupled dou- 
ble dot is placed in one arm of a "two-path" Aharonov- 
Bohm interferometer and the current oscillations as a 
function of a magnetic flux enclosed by the two arms are 
measured. Such measurements were already performed 
for single multi-level dotsi24 Across a transmission reso- 
nance the phase continuously changes by tt. Associated 
to G{Vg = 0) = is a phase jump of tt. From the 
conductance can be computed as Gjie^jh) ~ 2sin^ao-- 
Note that this relation together with Eq. (|39|l only im- 
plies aa = ((Hct) +m)'K, with m d Z and the observed 
phase jumps by tt are not excluded. The weak shoulder 
at ao- « 7r/2 for the parameters of Fig. ^| is related to 
the Kondo effect and develops into a plateau if U /F is 
increased. 

From the Ti. = conductance curves the two regimes 
(local spin-singlet and two-stage Kondo) cannot be dis- 
tinguished unambiguously. In particular, in both cases 
G is suppressed very strongly close to half-filling. The 
G{Vg) curves in the two regimes change quite differently 
in the presence of a magnetic field that lifts the spin- 
degeneracy of both levels: £ij = £2.1 = — 'H/2, 
ei^l — £2a = Tg + H/2. As for the single dot case 
studying the magnetic field dependence also enables us 
to define Kondo temperatures. For the two parameter 
sets of Figs. El and G{Vg) for different H is shown 
in Figs. 1131 and 1141 In the two-stage Kondo regime we 
define a scale T]^ for the first stage of the Kondo ef- 
fect considering the single dot case, i.e. setting ti^2 — 0. 
For the parameters of Fig. E] we obtain T^/F — 1.16 
in good agreement with the NRG estimate T^/T = 1.^- 
Because of its very small energy scale Tq <C the sec- 
ond stage of the Kondo effect gets destroyed already by 
a tiny magnetic field Ti. <C and the local minimum 
of G at half-filling evolves into a maximum as shown in 
Fig. We define Tq by the magnetic field required to 
obtain G{U,Vg = 0,H = To) = e^/h. For the param- 
eters of Fig. El this leads to Tq/F = 0.003. For larger 
ti,2 (spin-singlet regime), as shown in Fig. 1131 the spin- 
singlet is the ground state up to moderate Ti and Vg = 
remains a minimum. The magnetic field dependence of 




FIG. 13: (Color online) Gate Voltage dependence of the total 
conductance G for the parameters of Fig. IIUI and different 
magnetic fields H ^ (solid line), H = 0.065r = 0.5Tk 
(dashed line), Ti. = O.ISF = Tk (dashed-dotted line). 
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FIG. 14: (Color online) Gate Voltage dependence of the to- 
tal conductance G for the parameters of Fig. llll and different 



magnetic fields Ti = (solid line), Ti. 
0.0026r]^ (dashed line), H = O.OSF = 
(dashed-dotted hne). 



o.oosr 



ior° = o.026r]^ 



the two resonances is similar to the single dot case. In 
this parameter regime it is meaningful to define a Kondo 



temperature Tk by G{U, Vg 



±VgP,n 



Tk) 



For the parameters of Fig. El we obtain Tk/Y — 0.13, 
which is close to the NRG resuff T^/F = 0.156.2« In 
both parameter regimes the evolution of G{y„) with Ti 
is similar to the one with T discussed in Ref. l28l 

For the present setup we also studied the role of (a) 
asymmetric dot- lead couplings ^ F{*, (b) asymmetric 
local interactions Ul'i ^ C^J'i ' ^^^d (c) a nearest-neighbor 
interaction V^'^ , with j denoting the complement of j. 
In case (a) the overall height of G decreases with increas- 
ing asymmetry while the line shape remains invariant 
similar to the case of a single dot.^** Asymmetric local 
interactions and nearest-neighbor interactions only lead 
to quantitative changes but no new physics. 



VII. PARALLEL DOUBLE DOTS 

As our most complex example we finally study the par- 
allel double dot of Fig. ^E. The large number of param- 
eters allows for a variety of different regimes and we here 
refrain from giving a complete account of the physics of 
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this setup. Instead, we focus on a few cases to exemphfy 
that our approximation scheme can be used to uncover 
interesting effects. A more detailed discussion will be 
given in an upcoming publication. We here only con- 
sider the case without an inter-dot hopping. If present 
initially such a term can always be tuned away by a 
basis transformation on bonding and anti-bonding dot 
states. Note that in the noninteracting case the double 
dot can equally be viewed as a single two-level dot. In- 
cluding the interaction both cases might be distinguished 
by the choice of non-vanishing matrix elements (see be- 
low). Correlated multi-level dots have recently attracted 
much interested in connection with the puzzling obser- 
vation of "universal" phase lapses in a series of measure- 
ments of the transmission phase through quantum dots.^^ 
The conductance for the present setup with full spin po- 
larization and small level detuning was earlier studied us- 
ing the simplest fRG-based approximation scheme (con- 
sidering only the flow of the self-energy) and NRG. In the 
absence of spin-Kondo physics the inter-dot (inter-level) 
interaction leads to novel correlation-induced resonances 
located exponentially (in U) close to Vg = Om^ 

The noninteracting projected "effective Hamiltonian" 
of the present setup is given by 

/^.(^c.) = ( si,- rf -^7Sgn(a;) \ 

^ ' V sgn(w) e2,<T - sgn(w) y ^ ' 



with 

- E r!- ' 7 = V^rfit + e'* y^rfif . (4i) 

We furthermore define 

Here (p denotes a magnetic flux that goes through the 
ring geometry. If the system is viewed as a single two- 
level dot, (f> is restricted to and tt, on which we will 
focus in the following. We denote the relative sign of the 
hopping matrix elements i^- by s = e**^ = ±1. The non- 
interacting propagator Go(i^) follows by the matrix in- 
version Eq. H3()l) . For generic parameters during the fRG 
flow Eqs. IjlQI) and (|20() the level positions are renormal- 
ized and a (real) direct inter-dot hopping is generated. 
Applying the projection technique Eq. H1U|) the conduc- 
tance can be computed from the matrix elements of the 
projected Green function at A = 



h 



rf rf c;^,T°i(o) + V^^rf t;^-i(o) + Wrf r^^tj^^Ko) + Wr^r^^t;^-i(o) 



(43) 



The large number of parameters make it essential to analyze the conductance for the noninteracting case before 
considering the effect of two-particle interactions. A closed expression for Ga-{Vg) at Up^ — can be obtained from 
Eq. 113 by replacing ^^^="(0) by ^^0(0) 



GAVg) = -r 



4 










2 

+ [ei.<Tr2 + 





(44) 



We here focus on the spin-degenerate case with 7i = 0. 
For generic level-lead couphngs'^^ F^ and 5 = £2,0- — £1,0- > 
the gate voltage dependence (£1/2,0- — + Vg) of 

Eq. H44|) is characterized by two peaks (of height < /h) 
and a conductance zero. The latter follows from perfect 
destructive interference at a particular Vg. Associated 
to the zero is a jump of the transmission phase by tt. 
Further details of Ga{Vg) depend on s. For s = -1-1 the 
conductance zero is located between the two conductance 
peaks for all 5. For 5^0 the peak positions depend on 
the asymmetry of the F^- and the separation of the peaks 
is small if the F^- are almost equal. For equal F^- and (5 = 
the zero disappears. For s — —1 and fixed F^ the posi- 
tion of the conductance zero with respect to the peaks 



depends on 5. For 5 ^ it is located between the two 
conductance peaks, while it lies outside for large 5. In 
the crossover regime between these limiting cases one of 
the peaks vanishes, while the other becomes broader and 
splits up into two resonances separated by a minimum 
with non-vanishing conductance. For F^- with 1-2 sym- 
metry and (5 — > the conductance vanishes for all Vg as 
the model can be mapped onto a model with two levels, 
each one only coupled to one of the leads. In this (and 
only this) case a pseudo-spin variable is conserved leading 
to an orbital Kondo effect if a nearest-neighbor interac- 
tion is added. '^^■^^ For fixed, asymmetric F^- and 5 Q 
the separation of the two conductance peaks for s — —1 
is much larger than for s = +1. For certain classes of 
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FIG. 15: (Color online) Gate Voltage dependence of the total 
conductance G of a parallel double dot (see Fig. Q E) with 
generic Fj, spin-degenerate levels, purely local interactions, 
and different level detunings S. The parameters are U/T — 4, 
Ff/F = 0.27, Ff/F = 0.33, F^-ZF = 0.16, Ff/F = 0.24, 
s = 1, 5/F = (solid line), S/F = 1 (dashed line), and <5/F = 5 
(dashed-dotted line). 



parameters, some of them we mentioned already (e.g., 
equal and S = 0; see Ref. IssI for details), G{Vg) does 
not follow the behavior described above. Here, we do not 
investigate these special cases but focus on the generic 
behavior. 

In Fig. El we show G{Vg) for a generic set of F^ , 
s = +1, and different S. We assume a purely local inter- 
action U, which might experimentally be realized if the 
setup is implemented as two single-level dots of sufficient 
spatial separation. Close to half-filling the conductance 
is strongly suppressed. In particular this holds for small 
S. We note that around Vg = the flowing two-particle 
vertex becomes large which indicates that our approxi- 
mation scheme becomes less accurate. Further down we 
return to this issue. For S/T < 1 the curves of Fig.lTSlare 
similar to the one of Fig.llllobtained for the side-coupled 
double dot with small intcr-dot hopping ti^2- We argue 
that this similarity follows from a relation between the 
two double-dot geometries which we investigate next. 

Taking the even and odd linear combination of dot 
states the noninteracting side-coupled dot (index s) with 
(ti)^, (ti)^, {ti,2)s, and (Sj^a)^ = Vg can be mapped on 
the parallel dot with parameters 

ei.a ^ Vg- (ii,2), , e2.a ^ Vg + (ii,2), , (45) 

that is, 5 = 2(^1.2)3. Note that small (^1,2)3 correspond 
to small 5. This subset of the large parameter space of the 
parallel double dot is characterized by a 1-2 symmetry of 
the hoppings and s = 1. A local interaction of the side- 
coupled dot 

f^E-^.T"L- (46) 
3 



maps onto 

+ C?l,|rf2,T'^2,J,c;i:i + rfI,|'^2,T'^l,J.'^2,i) (47) 

for the parallel double dot. Besides the usual density 
interactions also correlated hoppings are generated by 
this transformation^ 

These considerations show that the cases studied in 
Figs. 111! and 1151 (for 5/T < 1) are related but not com- 
pletely equivalent as in Fig. El we only considered a lo- 
cal interaction and no 1-2 symmetry of the F^-. Slowly 
increasing the amplitudes of the additional interactions 
Eq. (|47|l and tuning the F' towards 1-2 symmetry we 
convinced ourselves that both the symmetry and the ad- 
ditional interactions do not play an essential role and the 
results for the parallel double dot with S/T < 1 of Fig. El 
can be interpreted as for the side-coupled double dot in 
the two-stage Kondo regime. The strong suppression of 
the conductance of the parallel double dot around half- 
filling for small 5 then follows from the second stage of 
the Kondo effect. From the side-coupled dot we already 
know that our approximation scheme overestimates this 
suppression and we believe that the same holds for the 
present geometry. Note that for (5 = 0, G(Vg) is symmet- 
ric around = as the Hamiltonian is symmetric under 
a particle-hole transformation if in addition Vg — )■ —Vg. 
Despite the asymmetry of the F^- the two resonances have 
almost unitary height. For small 6 the occupancies of the 
two levels are similar, sharply change across the two res- 
onances, and display a plateau at 1 close to = 0. The 
transmission phase increases by tt across the resonances 
and jumps by tt in between the peaks. 

For large S, G{Vg) in Fig. El can be understood from 
the eigenvalues of the many-body dot Hamiltonian with- 
out the leads, in analogy to the limit of large (^1.2)^ for 
the side-coupled dot. Each of the resonances then corre- 
sponds to a single-level resonance, as discussed in Sec. lIVI 
Because of the asymmetry of the F^-, the conductance 
does not reach the unitary limit and the two resonances 
have different line shape. Associated with the resonances 
is a change of one of the {rij) from 2 to (change of the 
phase from to tt) and a plateau at 1 (at 7r/2 for the 
phase) within the resonance. For the F^ of Fig. El ("-2) 
changes across the peak at < and (rii) across the 
one at Vg > 0. The plateaus in (rij) and are due to 
the Kondo effect that is active at each of the resonances. 
At the conductance zero between the conductance peaks 
the transmission phase jumps by tt. 

Scanning the parameter space with fixed s = -1-1 we 
encountered an interesting correlation effect for a model 
with almost equal local and nearest-neighbor density in- 
teractions, small 6, and generic F' . It is exemplified in 
Fig. El In the lower panel of Fig. 1161 we show G{Vg) for 
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FIG. 16: (Color online) Gate Voltage dependence of the to- 
tal conductance G of a parallel double dot with generic F' , 



spin-degenerate levels, and (almost) equal local and nearest- 
neighbor density interactions. The parameters are Ff/F = 
0.5, Ff/F = 0.25, F^/F = 0.07, /T = 0.18, and s = 1. 
Lower panel: Equal level positions on both dots S — 
and equal local and nearest-neighbor interactions U with 
U/T = 0.5, (solid line) 1 (long dashed line), 2 (dashed-dotted 
line), and 4 (short dashed line). Upper left panel: Stabil- 
ity of the novel resonances for non-vanishing level detuning 
S at U/r = 2. 5/T = 0.03 (solid line) and 5/V = 0.07 
(dashed line). Upper right panel: Stability of the novel reso- 
nances for asymmetric interactions at 5 = 0. U^'^ /T = 1.9, 
Ul'i/T = 2.1, Ujf/F = 1.8, Ujf/T = 2 (solid line) and 
U^f/T = 1.8, Ulf/F = 2.2, C/ff /F = 1.7, U^f/T = 2.2 
(dashed line). Note that here 



3,3 ' 3,3 

is shown in units of F. 



i5 = and U"'^ = U{1 — Sjj'S„^„'). Increasing U, the 
peak conductance increases, the two maxima become flat, 
and for U larger than a F^- dependent critical Uc{{T'j}) 
both split up into two peaks. For U > Uc the height 
of all four resonances is equal and U independent. The 
two outermost peaks are located at « ±C/ and the 
inner two at ~ ±/7 /2. This scenario is quite similar 
to the appearance of additional correlation induced reso- 
nances in a spin-polarized model of a parallel double dot 
with nearest-neighbor interaction as discussed in Ref. IsM 
The novel resonances follow from the combined presence 
of correlations and quantum interference. In the upper 
panels of Fig. E| we show that the effect is robust against 
a small level detuning d and a small asymmetry of the 
JJjy . It will be further investigated in an upcoming pub- 
lication. One can expect that in experimentally realized 
single two-level dots the inter- and intra-dot density in- 
teractions are indeed (almost) equal and the additional 
resonances should thus be observable in such systems. 

We next study the case s = — 1. As for s = -f 1 above 
we first investigate systems with purely local interactions 
(two well-separated single-level dots) and further down a 
situation with (almost) equal local and nearest-neighbor 
interactions (a single two- level dot). 




FIG. 17: (Color online) Gate Voltage dependence of the to- 
tal conductance G of a parallel double dot with generic Fj, 
spin-degenerate levels, purely local interactions, and different 
level detunings 5. The parameters are U/T — 4 (thick lines), 
U/F = 1 (thin line) Ff/F = 0.5, Ff/F = 0.35, F^/F = 0.05, 
F^/F = 0.1, s = -1, 5/F = (solid line), 5/F = 0.6 (dashed 
line), and 5/F — 6 (dashed-dotted line). 



In Fig. El we show G{Vg) for s = —1, generic F^ , a 
purely local interaction, and different 6. Without loss of 
generality we consider F^- with a fairly large 1-2 asymme- 
try. Otherwise the conductance at small S is very small. 
This follows from the fact that for 1-2 symmetric F^- and 

5 — 0, G{Vg) — for all Vg, as discussed above. The thin 
solid curve shows the case with 6 = and U/T = 1, that 
is qualitatively similar to the U = case. Increasing U 
(at 6 = 0) each of the two peaks splits up into two with 
a conductance zero in between. In Fig. 1171 this is shown 
as the thick solid line for U/T = 4. The dashed and 
dashed-dotted lines show how the conductance evolves 
with increasing level detuning 6 at fixed U/T = 4. The 
conductance zeros either vanish or are pushed towards 
larger \Vg\. The overall conductance increases. For large 

6 the two remaining resonances can be understood as 
in the single-dot, single-level case. Across the resonance 
centered around —6/2 (+6/2) level 2 (1) gets filled. As 
the F2 are small the effective interaction U/{T2 + T2) 
relevant for the Kondo effect is large and the resonance 
at < has a well-developed plateau-like line shape. 
These results can most easily be understood by making 
contact with the case of left-right symmetric F' . In this 
case the off-diagonal element 7 [see Eq. (|41|) ] of the non- 
interacting "effective Hamiltonian" 140() vanishes. As the 
interaction is purely local no such term (that is a hop- 
ping between j = 1 and j = 2) is generated during the 
fRG flow and the double-dot problem separates into two 
single-dot problems. In particular, Ga^% vanishes and 
the conductance Eq. H43|) is given by the sum of two 
single-level Green functions. Breaking of the left-right 
symmetry leads to a small initial off-diagonal element, 
that stays small during the fRG flow. For that reason also 
the results of Fig. El can be understood from the single- 
level case. E.g., this explains the position Vg ~ ±6/2 of 
the resonances at large 6. A more detailed discussion will 
be given elsewhere. 

In Fig. 1181 we show how the four-peak structure ob- 
tained for a sufficiently large local interaction U at 6 = 
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FIG. 18: (Color online) Gate Voltage dependence of the total 
conductance G of a parallel double dot with generic F^ , s = 
— 1, spin-degenerate levels, 5 = 0, local interaction U/T = 4, 
and different nearest- neighbor interactions [/' = (solid line) , 
U'/r = 2 (dashed line), and U' /F = 4 (dashed-dotted line). 
The Fj are as in Fig. 1171 Note that here Vg is shown in units 
of F. 

and s = — 1 disappears in the presence of a nearest- 
neighbor interaction U' (single two- level dot). Starting 
out with four peaks of almost equal height at U' — 0, 
for increasing U' the two outermost peaks gain weight, 
while the inner two peaks lose weight. At U' = U only 
two resonances remain. 



VIII. SUMMARY AND OUTLOOK 

In the present paper we have introduced a novel ap- 
proximate method to study electronic transport through 
systems of quantum dots with local Coulomb correla- 
tions. Our fRG-based approximation scheme results 
in a set of coupled differential flow equations for the 
interaction-dependent effective level positions and inter- 
level hoppings as well as the interaction matrix elements. 
We have solved the flow equations for a variety of dot 
systems and computed the linear conductance, the trans- 
mission phase, and level occupancies. 

For several setups and parameter sets we compared 
our results to existing exact Bethe ansatz and high pre- 
cision NRG data. For a single dot the agreement is ex- 
cellent up to the largest Coulomb interaction for which 
Bethe ansatz or NRG data are available. In particular, 
we showed analytically that our approach covers the sig- 
natures of the Kondo effect. The method performs very 
well also in the presence of a magnetic field lifting the spin 
degeneracy and in case of an additional direct transmis- 
sion channel between the leads (Fano effect). The field 
dependence of G{Vg) allows us to extract the Kondo tem- 
perature that compares well to the exact expression. For 
linear chains of dots (short Hubbard chains) we found 
quantitative agreement with NRG data up to fairly large 



interactions. We studied the effect of a left-right asym- 
metry of the level-lead couplings that leads to the sur- 
prising vanishing of resonances if sufficiently large. For 
a side-coupled double-dot system we demonstrated that 
our method can be used to describe both the "molecu- 
lar" spin-singlet regime at large inter-dot hoppings ti^2 
and the two-stage Kondo regime at small ti,2- For not 
too small ti^2 and too large interactions the approximate 
results agree well with NRG data. Our approximation 
scheme becomes questionable if C//i^ 2 becomes too large. 
We showed that the occupancies in the two-stage Kondo 
regime depend non-monotonically of the gate voltage. At 
T — the different physics of the two regimes can be 
probed by a magnetic field lifting the spin degeneracy. 
Finally we demonstrated that the physics of a parallel 
single-level double dot (single two-level dot) is very rich. 
For a certain parameter regime we showed that G{Vg) can 
be understood in analogy to the side-coupled double dot. 
For other classes of parameters we surprisingly found ad- 
ditional resonance peaks that should be observable exper- 
imentally in existing double-dot setups (two- level dots). 
A more complete analysis of the transport properties of 
the parallel double dot (single two-level dot) will be given 
in an upcoming publication. 

The present paper shows that our approximation 
scheme is a simple-to-implement and reliable method to 
study transport through dot systems with a significant 
number of interacting degrees of freedom. It can be used 
to study setups that due to their complexity cannot be 
tackled by NRG. Even in problems that can be investi- 
gated by NRG our method is far superior in terms of the 
numerical effort required, e.g. enabling efficient analysis 
of parameter dependencies. 

In summary, we have shown that the fRG based ap- 
proximation scheme is a very promising tool to inves- 
tigate the interesting physics of quantum dot systems 
resulting from Coulomb correlations (Kondo effect) and 
quantum interference. 
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